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Abstract 

In the scenario of isospin- violating dark matter (IVDM), the dark matter (DM) 
spin-independent couplings to protons and to neutrons are allowed to be different, 
which has been considered to relax the tensions between the results of DAMA, 
CoGeNT and XENON experiments. We explore the allowed values of DM-nucleon 
couplings favored and excluded by the current experiments under the assumption 
of IVDM. We find that the recently updated XENON100 result excludes the main 
part of the overlapping signal region between DAMA and CoGeNT. The mass of the 
IVDM is restricted to a narrow range near 7.5-8 GeV. We also show that the possi- 
ble tensions between some experiments such as that between DAMA and SIMPLE 
are unlikely to be affected by isospin violating interactions. In an effective operator 
approach, we investigate conservative upper bounds on the DM-quark couplings 
required by the IVDM scenario from the cosmic ray antiproton fluxes measured 
recently by BESS-Polar II and PAMELA, and that from the relic density. The 
results show that the relatively large couplings favored by DAMA and CoGeNT 
are tightly constrained, if the operators contribute to velocity-independent anni- 
hilation cross sections. For thermal relic DM, the upper bounds from the relic 
density can also be stringent. 
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1 Introduction 



It has been well established from astrophysical and cosmological observations that nearly 
85% of the matter in the Universe consists of invisible dark matter (DM). So far the 
evidence of DM arise solely from its gravitational interactions. Popular DM candidates, 
such as the weakly interacting massive particles (WIMPs) can naturally reproduce the 
observed relic density through weak interactions with the standard model (SM) particles. 
This possibility has motivated numerous experiments to probe the direct or indirect 
signals of DM interacting with ordinary matter. 

Some of the recent DM direct detection experiments such as DAMA [TH3], Co- 
GeNT jUE] and CRESST-II [6] have reported events in excess of known backgrounds. 
The excess events, if interpreted in terms of DM particle elastic scattering off target nu- 
clei, may imply light DM particles with mass around 8-10 GeV and scattering cross sec- 
tion around 10~ 40 cm 2 . Other experiments such as CDMS-II [318], XENON10/100 0H}], 
and SIMPLE [11] etc., have reported null results in the same DM mass range. A com- 
monly adopted assumption on interpreting the DM direct detection data is that in spin- 
independent scatterings the DM particle couplings to proton (f p ) and to neutron (/„) 
are nearly the same, i.e. f n ~ f p , which makes it straight forward to extract the DM- 
nucleon scattering cross sections. However, in generic cases, the interactions may be 
isospin- violating [T2TT17] . In the scenario of isospin violating DM (IVDM), the DM 
particle couples to proton and neutron with different strengths, possible destructive in- 
terference between the two couplings can weaken the bounds from XENON10/100 and 
move the signal regions of DAMA and CoGeNT to be closer to each other JT6JEL7]- In 
order to reconcile the data of DAMA, CoGeNT and XENON10/100, a large destructive 
interference corresponding to f n / f p ~ —0.7 is favored [T6] . 

Recently, the XENON100 collaboration has reported updated results based on the 
exposure of 225 days with 34 kg fiducial mass [18]. Two events are observed with 
background expectation of 1.0 ± 0.2 events (0.79 ± 0.16 from gammas and 0.17+°; J 2 from 
neutrons). The updated upper bound on the spin-independent scattering cross section is 
2 x 10 -45 cm 2 at 90% confidence level (CL) for a 55 GeV DM particle, which is improved 
by a factor of 3.5 in comparison with the previous result reported in 2011 [10]. For 
an 8 GeV DM particle, the limit is improved by a factor of five. Such a significant 
improvement may allow the XENON100 result to challenge the DAMA and GoGeNT 
results even in the scenario of IVDM. 

As a consequence of f n / f p ~ —0.7, the absolute values \f p , n \ have to be around two 
order of magnitudes larger than that obtained under the assumption of IC interaction. 
At quark level, this means that the absolute DM couplings to the first generation u- and 
G?-quarks must be enhanced by the same order of magnitudes. The scenario of IVDM 
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with such light DM masses and relatively large couplings to light quarks can be tested in 
other processes in which there is no destructive interference. For instance, the light DM 
particles can be pair produced at the Large Hadron Collider (LHC) and the Tevatron 
with signals of a single jet/photon plus a large missing transverse energy, i.e. a mono- 
jet/mono-photon. The current null search results from LHC [T9l[20] and Tevatron [21] 
can be used to impose constraints on the DM-quark couplings [22TT25] and the IVDM [26] . 

As the inverse processes of DM production, light DM particles can annihilate into 
light-quark pairs in the galactic halo, which provides exotic sources of cosmic gamma 
rays, neutrinos and antiprotons. The annihilation processes occur at low velocities v/c ~ 
O(10~ 3 ) where c is the speed of light, which can be used to test the IVDM complementary 
to that at the LHC and Tevatron. For such low energy processes, the effective operator 
approach can be valid as long as the intermediate mediator particles are significantly 
heavier than the DM particle, but could be much lighter than the TeV scale (for recent 
discussions, see e.g. Refs. [23 | [2 ^|27H30] ) . Possible constraints from the cosmic neutrinos 
and gamma ray on IVDM have been discussed in Refs. [3TH33] . Recently the BESS-Polar 
II experiment has measured the antiproton flux in the energy range from 0.2 GeV to 3.5 
GeV [M] which have higher precision compared with that from PAMELA [35] at low 
energies. The antiproton flux in this energy range can receive significant contributions 
from the annihilation of 10 GeV scale DM particles, provided that the annihilation cross 
section is not velocity suppressed. The constraints from BESS-Polar II data have been 
investigated recently for typical thermal WIMP [36] . In this work, we shall focus on the 
implications of the recent cosmic ray antiproton flux data for the scenario of IVDM and 
use different approach to calculate the propagation of anitproton. 

For thermal relic DM, the same annihilation process also determines the DM relic 
density. By requiring that the calculated relic density from a single annihilation channel 
should not be smaller than the observed value flh 2 = 0.113 ± 0.004 [37], upper bounds 
on the DM couplings can be obtained. The bounds could be stringent if the annihilation 
is dominated by s-wave processes. For p-wave annihilation, although the cross section 
is velocity-suppressed, useful upper bounds can still be obtained as the typical relative 
velocity of DM particles is finite v jc m 0.3 at freeze out, 

In this work, We first explore the values of DM-nucleon couplings favored by the 
current experiment under the assumption of IVDM for different target material and 
values of f n /f p - We find that the recently updated XENON100 result is able to rule 
out the main part of the overlapping signal region between DAMA and CoGeNT. The 
mass of the IVDM is restricted to a narrow range near 7.5-8 GeV. We also show that 
the tensions between some group of experiments are unlikely to be affected by isospin 
violation, especially that between DAMA and SIMPLE. We adopt the effective operator 
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approach to investigate the constraints from the antiproton flux data and the thermal 
relic density on the couplings between the IVDM and the SM light quarks. We calculate 
the antiproton flux using the numerical method implemented in the GALPROP code [381 - 
|4"T] . and consider a number of propagation parameter configurations in order to obtain 
conservative upper bounds. The results show that the large and negative value /„/ f p = 
—0.7 is severely constrained for the operators with velocity- independent annihilation 
cross sections, such as the fermionic DM with vector couplings and complex scalar DM 
with scalar couplings. For these operators we also find that the constraints from the 
cosmic antiproton flux are stronger than that from the thermal relic density. The relic 
density can provide useful constraint in the case where the operators contribute only 
to velocity suppressed annihilation cross sections. For the complex scalar DM with 
derivative couplings, we find that the IVDM is in some tension with the relic density. 

This paper is organized as follows. In Section |2} we discuss the effect of IVDM on 
various target nuclei in different experiments and then explore the allowed DM couplings 
to nucleons and to light quarks. In Section |3l we list the effective operators relevant to 
the IVDM and the related interaction cross sections. In Section HJ we derive constraints 
on IVDM from the current cosmic ray antiproton data. In section |5] the constraints from 
the DM relic density is discussed. The numerical results are presented in Section [6] and 
the conclusions are given in Section [7J 



2 Isospin violating dark matter 

For a DM particle \ with mass m x elastically scattering off a target nucleus with atomic 
number Z and atomic mass number A, the recoil event rate R = dN/dt is given by 

where Nt is the total number of the target nuclei and po is the local DM energy density. 
For a fixed recoil energy Er, the required minimal velocity of the DM particle is v m i n = 
[itiaEr/ (2/i^)] 1//2 , where vtla is the mass of the target nucleus and fx a = ( m x m A> / ( m x + 
ttia) is the DM-nucleus reduced mass. The maximal velocity is the escape velocity from 
our Galaxy at the position of the Solar system v esc ~ 600 km ■ s" 1 . The differential 
scattering cross section can be written as 

da m A E 2 (E R ) 



dE R 2fx 2 A v 2 

where F(Er) is the form factor of the nucleon. Here we have assumed that the form 
factors for proton and neutron are nearly identical, i.e., F p {E R ) « F n (En) = F(Er). 
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The quantity o"o can be understood as the total scattering cross section at the limit of 
zero-momentum transfer which is related to f p ( n \ through 

a = ^[Zf p + (A-Z)f n ] 2 . (3) 

7T 

The cross section for the DM particle scattering off a free nucleon in term of f p < n ) is 

where fi p ( n ) is the DM-proton (neutron) reduced mass. Under the assumption that the 
scattering is isospin conserving (IC), i.e., f n ~ f p , the total cross section <7o is independent 
of Z and only proportional to A 2 . One can define a cross section <j IC which is the value 
of a p extracted from a under the assumption of IC interaction as 

which is the quantity commonly reported by the experiments. In the generic case where 
fn 7^ f P , the true value of a p will differ from a p c by a factor F(f n / f p ) which depends on 
the ratio f n j f p and the target material 

op = F{f n lf p )af. (6) 

Depending on the mass of the DM particle, for a given target material, for instance 
CaW04 used by CRESST-II experiment there could be more than one kind of nuclei 
relevant to the nuclear recoil. If the target material consists of N kind of relevant nuclei 
with atomic numbers Z a [a = 1, . . . , N) and fractional number abundances n a , and for 
each nucleus Z a there exists M type of isotopes found in nature with atomic mass number 
A ai and fractional number abundance r] ai (i = 1, . . . , M), the expression of F(f n / f p ) can 
be explicitly written as 

r {Jn/Jp) = ^ - ,, 2 rv , , A ry —f 7Tp ' 



E tt ,i K aVai^ 2 A ai [Z a + (A ai Z a )f n / f 



where fiA ai is the reduced mass for the DM and the nucleus with atomic mass number 
A ai . In the simplest case where the target consists of one kind of nucleus (Z, A), F(f n / f p ) 
takes the simple form 

-2 



F(fn/f P ) 



A \ AJ f p 



(8) 



It is evident that F(f n / f p ) approaches unity in the case of IC scattering. If f n j f p <0 the 
interference between the contributions from proton and neutron scattering to the value 
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of F(f n / fp) will be destructive, which can lead to F(f n / f p ) ^> 1. Thus it is possible 
that the value of a p can be a few order of magnitudes larger than <J p c , provided a nearly 
complete cancellation between the two contributions. 

For a given target material T, there is a particular value of f n / f P which corresponds 
to the maximal possible value of F(f n /f p ) 

r Ea,! K aVai^A ai (Aai ~ Z a )Z a 

= 2 T~A V \2 ' ^ ' 

For a single nucleus target with atomic (mass) number Z(A), it is simply given by 
£z = —Z/(A — Z). The value of £t varies with target material. In Tab. [H we list the 
values of £t for some typical material utilized by the current or future experiments. 

Xe Ge Si Na(I) Ca(W)0 4 C 2 C1F 5 Csl Ar 
£r -0.70 -0.79 -1.0 -0.92(-0.73) -1.0(-0.69) ^92 -0.71 -0.82 

Table 1: Values of £t for different target material. For Nal, the two values -0.92 and 
-0.73 correspond to the scattering off Na and Nal respectively. Similarly, for CaWO"4, 
the two values -1.0 and -0.69 corresponds to the scattering without and with tungsten 
respectively. 



In Fig. [H we plot the regions favored and excluded by the current experiments 
in the plane of (f p , f n / f p ) for four different DM masses from 7.5 GeV to 9.0 GeV. 
The experimental results include: CDMS-II 90% CL upper bounds pJIB], CRESST- 
II 2a favored region for scattering without involving tungsten nucleus [6], DAMA 3a 
favored region assuming Na scattering without including the channeling effects [3], 90 
% CL upper bound from SIMPLE [IT], the 90% CL upper bound from XENON10 
which is insensitive to scintillation efficiency [9] and the 90% CL upper bound from the 
XENON100 p]. As it has been noticed previously at f n /f p = -0.70 = £ Xe , the DAMA- 
and CoGeNT-favored regions overlap, the sensitivities of XENON10/100 are maximally 
reduced by two order of magnitudes. However, the DAMA- and CoGeNT favored region 
is only marginally consistent with the new XENON limit for m x = 7.5 and 8.0 GeV. For 
m x = 8.5 and 9.0 GeV, there is no overlapping region. In Fig. [21 The allowed regions 
in (a p , m x ) plane are shown. One can see clearly that the main part of the overlapping 
region between DAMA and CoGeNT is excluded by the XENON100 result. 

If the £r values of the target material used by two experiments are very close to 
each other, the tension between the two experimental results, if exists, is less affected 
by the effect of isospin violation. From Tab. [U one finds that <^Na ~ £c 2 cif 5 = —0.92, 
6ce ~ £csi ~ —0.7 and £si ~ £ca(w)o 4 = —1-0. Thus the tension between DAMA signal 
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from Na recoil and the upper bound from SIMPLE is unlikely to be alleviated by isospin 
violation, which can be clearly seen in Fig. [TJ Similarly, if there exists contradictions 
between XENON and KIMS, CoGeNT and the Ar based experiments such as DarkSide, it 
can hardly be explained by isospin violating scattering. The SIMPLE result is also useful 
in comparing with the CRESST-II which utilizes Ca(W)0 4 which has £ca(w)o 4 — —1-0. 
Obviously, for the experiments use the same target material, the possible tension between 
them cannot be relaxed by isospin violation, such as the tension between CoGeNT and 
CDMS-II, as both use germanium as target nucleus. 
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Figure 1: Favored regions and upper bounds in (f p , f n / f P ) plane from current experi- 
ments DAMA [2], GoGeNT [5], XENON PUB], CDMS (HE] and SIMPLE [II] are also 
shown. Four panels corresponds to the four different mass of dark matter particle fixed 
at 7.5 GeV, 8.0 GeV, 8.5 GeV and 9.0 GeV respectively. 



3 Effective interactions 

We assume that the DM particles interact with the SM light quarks through some heavy 
mediator particles much heavier than the DM particle such that both the scattering 
and the annihilation processes can be effectively described by a set of high dimensional 
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Figure 2: The favored regions and constraints in the (a p , m x ) plane for various experi- 
ments for f n / f p = —0.70. 

contact operators 

£ = Y,aiA- (10) 

i,q 

If the DM particles are Dirac fermions, the relevant operators arising from scalar or 
pseudoscalar interactions are given by 

Oi g = XXQQ, 2q = X7 5 X<?<?, 3q = xXQl 5 (l, O iq = xi^XQl^Q- (H) 
The operators from vector or axial-vector type interactions are 

o 5q = xi^xti^ o &q = xYi 5 xqj»q, o 7q = xYxqjrfq, o 8q = xtVxw 5 ?, (12) 

and the ones from the tensor interactions are 

o 9q = x^xqw, £>m = x^^xqv^q- (is) 

If the DM particles are Majorana particles the vector and tensor operators are vanish- 
ing. Among these operators only G lq and 5q contribute to spin-independent scattering 
cross sections at low velocities. The scattering cross sections induced by the operators 
Oi q and 0<a q are velocity suppressed. The operators 7q and Os q contribute only to 
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spin-dependent scattering cross section, and the nucleus matrix elements for the oper- 
ators Ozq, O^g, Og q and Oiog are either vanishing or negligible. The DM annihilation 
into quarks through 0\ q is a p-wave process, which is velocity suppressed. It does not 
contribute to the cosmic antiproton flux, but still contributes to the DM relic density as 
p-wave processes is non-negligible at freeze out. 

Similarly, for DM being a complex scalar </>, possible operators are 

O n = 2m^(0*0)gg, 12 = 2m </ ,(0*0)g 7 5 g, 13 = (</>*% </>)<?fq, O u = (0*^ <f>)qj^ 5 q. 

(14) 

Among those only On and 13 contribute to the spin-independent scatterings. The DM 
annihilations through operator (9 13 are p-wave processes. 

The fermionic or scalar DM may couple to gluons through the operators XxG^G^ ', 
X r y 5 xG f iisG' Jll/ , xxG^ v G^ v and xi^xG^G^ . We do not consider these operators as they 
do not contribute to isospin violating scatterings. In summary, we only consider the four 
operators 

Oiq, O bq , llq , and O r3q , 

which are relevant to IVDM. 

The DM couplings to nucleons f p>n can be expressed in terms of the DM couplings 
to quarks follows 

For the Dirac DM with scalar interaction aiqXXQQi one nas B^^ = Tn p ( n )/m q 
for q = u,d,s and B^ q = (2/27) f^c m p{n) / m q f° r Q — c ,b,t, where f^ 1 is the DM 
coupling to light quarks obtained from the cr-term (N\m q qq\N) = f^ q M N , and f^Q = 1 — 
T, q=u ,d,s fTq ] - In numerical calculations we take /£ u = 0.020±0.004, f Td = 0.026±0.005, 
f Ts = 0.118 ± 0.062, f» u = 0.014 ± 0.003, f$ d = 0.036 ± 0.008 and f$ B = 0.118 ± 
0.062 [42]. Using the following quark masses: = 0.005 GeV, m u /rrid = 0.55, m s = 
0.095 GeV,m c = 1.25 GeV, m b = 4.2 GeV and m t = 172.3 GeV, we obtain the values of 

Bi q . 

UP ~ H» r^j P. 8 RP ~ R 11 ~ A 8 
D \u — D ld — u -°> D \d — D \u — ^-°) 

B*? ~ 1.2, B p £ ~ 0.05, B p { b n ~ 1.5 x 10~ 2 , and ~ 3.5 x 10" 4 . (16) 

In order to maximize the isospin violating effect, the coefficients B\f lclblt must be 
strongly suppressed. Assuming that the DM-nucleon couplings are dominated by the 
DM couplings to the first generation quarks, the ratio f n / f p is given by 

fn ^ B r { u a lu + B'1 d a ld 

~~F~ ~ "dp dp ' v / 
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The value of f n / f P = —0.7 can be translated into a\dja\ u = —0.93 at quark level. This 
value is the same for complex scalar DM. For operator 0§ q one simply has f?^ = 2(1) 
and Blf ] = 1(2), and B p q [n) = for q = c, s, t, b. 

fn «5tt + 2<25d 



(18) 

Jp 2a5 U + a^d 

Thus for f n / f p = —0.7, one finds a^d/^bu = —0.89. The cross section for DM annihilating 
into quarks induced by the relevant operators are given by 



Ncm 



2 

VlqV le l(xX -> ft) = Z X a iq V reh 



2tt 

Ncm 



2 

x Jl 



71 



5g i 



O-llg^rcl(00 -> 55) = 



,2 

- \ 2iV ^ 2 2 / 1n \ 
^13^rcl(00 -> W = g QlSg^reb ( 19 ) 

where Nc = 3 is the number of color and v re \ is the relative velocity of annihilating 
DM particles. For s-wave annihilation, the thermally averaged value of the cross section 
multiplied by velocity (av Te \) is the same as av Te \. 



4 Cosmic antiproton flux 

Annihilation or decay of light DM particles in the galactic halo can contribute to ex- 
otic primary sources of the low energy cosmic ray antiprotons, which can be probed or 
constrained by the current satellite- and balloon-borne experiments such as PAMELA 
and BESS-polar II, etc.. The predicted antiproton flux from DM annihilation depends 
on models of the cosmic-ray transportation, the distribution of Galactic gas, radiation 
field and magnetic field, etc.. It also depends on the particle and nuclear interaction 
cross sections. In this work we use the numerical code GALPROP [38TI¥T[H3] which uti- 
lizes realistic astronomical information on the distribution of interstellar gas and other 
data as input and consider various kinds of data including primary and secondary nu- 
clei, electrons and positrons, 7-rays, synchrotron and radiation etc. in a self-consistent 
way. Other approaches based on simplified assumptions on the Galactic gas distribution 
which allows for fast analytic solutions can be found in Refs. [HH17] . In the GALPROP 
approach, we consider several diffusion models (parameter configurations). The different 
results between the models can be regarded as an estimate of theoretical uncertainties. 

In the diffusion models of cosmic ray propagation, the Galactic halo where diffusion 
occurs is parameterized by a cylinder with half height Zh and radius R = 20 — 30 kpc. 
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The densities of cosmic ray particles are vanishing at the boundary of the halo. The 
processes of energy losses, reacceleration and annihilation take place in the Galactic 
disc. The source terms for the secondary cosmic rays are also confined within the disc. 
The diffusion equation for the cosmic ray particle is given by 

-—i>--i/} + q(r,p), (20) 
Tf r r 

where ip(r,p,t) is the number density per unit of total particle momentum which is 
related to the phase space density f(r,p,t) as ip(r,p,t) = 4np 2 f(r,p,t). For steady- 
state diffusion, it is assumed that dip/dt = 0. The spatial diffusion coefficient D xx is 
parameterized as 

D xx = (3D (^j , (21) 

where p = pj (Ze) is the rigidity of cosmic ray particle and 5 is the power spectral index 

which may take different values 5± or 62 when p is below or above the reference rigidity 

Po- D is a normalization constant, and = v/cis the velocity of the cosmic ray particle. 

The values of D and 5 are determined from the measurements of the fluxes of other 

cosmic ray species such as the ratio of Boron to Carbon (B/C) and of isotopes 10 Be/ 9 Be. 

The convection term is related to the drift of antiproton from the Galactic disc due to the 

Galactic wind. The direction of the wind is usually assumed to be along the z-direction 

which is perpendicular to the galactic disc and is a constant Vc = \29(z) — 1]V C . The 

diffusion in momentum space is described by the reacceleration parameter D pp which is 

related to the Alfven speed V a of the disturbances in the hydrodynamical plasma as 

41/2 2 

n - VaP (02) 

Uw 3^(4-^)(4-^' [ZZ) 

where w stands for the level of turbulence which can be taken as w = 1 as D pp depends 
only on the combination V^/w. p is related to the momentum loss rate, 17 and r r are 
the time scale for fragmentation and radioactive decay respectively. 

The sources of the primary particles are chosen to reproduce the cosmic-ray distri- 
bution determined by the 7-ray data. The injection spectrum of nuclei is assumed to 
have a broken power low behavior dq(p)/dp ~ p~ 7 , with 7 = 71(72) for the rigidity p 
below (above) a reference rigidity p s . The secondary source is given in terms of the dis- 
tributions of the primary particles ip(r,p) and the distribution of the interstellar medium 
(ISM) 

q(r,p) = /3cVyi m (r,p)[a H (p)n H (r) + cr He nHe(r)], (23) 
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where nn(r) and %e(r) are the number densities of interstellar hydrogen and helium 
respectively, cth and cth g are the cross sections for the generation of the secondary 
particles from the interactions with H and He. The detailed calculations of the cross 
sections and the distribution of interstellar gas can be found in Ref. [39] and references 
therein. 

Thus the whole diffusion process depends on a number of parameters: R, Z h , p , D , 
81/82, V c , V a , p s and 71/72- There exists degeneracies in the determination of the diffu- 
sion parameters. In general, increasing the diffusion half-height Zh can be compensated 
by the increase of the diffusion constant Dq, which leads to significant uncertainties in 
the predictions of antiproton flux from DM annihilation. In some previous analysis ( see, 
e.g. Ref. [16]), the values of Z h can be as small as 1 kpc. The recent updated analysis 
show however that Zh should be around 4 kpc [4"7"1I4"8] . 

We consider three typical propagation models in GALPROP. The first one is the plain 
diffusion model (referred to as "Plain" ) in which there is no reacceleration term [JT] . The 
second one is the conventional model (referred to as "Conventional" ) with reacceleration 
included [401I4T]. The last one (referred to as "Global-Fit") is the model from a global fit 
to the relevant astrophysical observables using Markov-Chain Monte Carlo method [48J. 
The main parameters of the three models are listed in Tab. [2j Note that in the three 
models the convection is not considered, as the inclusion of conversion term leads to 
problems in reproducing the spectrum of B/C and discontinuity in the propagation 
across the Galactic disc in GALPROP [4T| . 



model 


^(kpc) 


Z^(kpc) 


£0 


Po 


61/62 


K(km/s) 


Ps 


71/72 


Plain 


30 


4.0 


2.2 


3 


0/0.60 





40 


2.30/2.15 


Conventional 


20 


4.0 


5.75 


4 


0.34/0.34 


36 


9 


1.82/2.36 


Global-Fit 


20 


3.9 


6.59 


4 


0.3/0.3 


39.2 


10 


1.91/2.40 



Table 2: Propagation parameters in the "Plain" [H] , "Conventional" [IQJHT] and 
"Global-Fit" 08] models used in the GALPROP code. D is in units of 10 28 cm 2 • s~\ 
the break rigidities p and p s are in units of GV. 

The primary source term from the DM annihilation has the form 

q(r) =r/n(r)V^ci>^, (24) 

where n(r) = p(r)/m x , 77 = 1/2(1/4) if the DM particle is (not) its own antiparticle and 
dN/dp is the injection spectrum per DM annihilation. For the DM profile we took the 
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isothermal profile [H] 

where r = 8.5 kpc is the distance of Solar system from the galactic center, R s = 2.8 kpc 
and the local density is taken to be po = 0.43 GeV cm -3 [50]. Note the predicted 
anitproton flux is modestly dependent on the halo profile. 

The total antiproton flux is related to the density function as 

$ = ^(p). (26) 

In the force- field approximation [51] , the antiproton flux at the top of the atmosphere 
of the Earth ^ OA which is measured by the experiments is related to the interstellar 
antiproton flux as 

% (T TOA )-{ 2mpT + T 2 )*p(n (27) 

where Ttoa — T — ipp is the antiproton kinetic energy at the top of the atmosphere of 
the Earth. The BESS-Polar II data were taken in a period of lowest solar activity. We 
take 4>p = 0.5 GV in numerical analysis. 



5 DM thermal relic density 

In the case where the DM particles are thermal relics, possible large couplings to light 
quarks may under predict the DM relic abundance in comparison with the observed 
value Vth 2 = 0.113 ± 0.004 [37]. Thus the relic density can also impose upper bounds on 
the relevant couplings. Whereas the p-wave annihilation processes give no contribution 
to the antiproton flux, in the calculation of DM relic density, the p-wave annihilation 
is nonnegligible as the typical relative velocity of DM particles is v Te \/c ~ (9(0.3) at 
freeze out. The annihilation cross section times the relative velocity can be expressed 
as av rd = a + bv^ e i, where a and b are coefficients corresponding to the s-wave and p- 
wave contributions. The thermally averaged value at the time of freeze out has the form 
(c^rei)/ — a + 6b/xj, where Xf = m DM /Tf with Tf the decoupling temperature. The 
value of Xf can be estimated through the iterative solution of the equation [52] 



Xf = In 



r>(r><o\ 9 M pi m DM(a + 66/37) 



(28) 



where g is the number of degrees of freedom of the DM particle, g* is the number of 
effective relativistic degrees of the thermal plasma and M p i = 1.22 x 10 19 GeV is the 
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Planck energy scale. The constant C is determined by matching the semi-analytical 
solutions to the fully numerical solutions to the Boltzmann equation for the thermal 
evolution of particle number density. We take C(C + 2) = 1(2) for s(p)-wave annihilation 
in numerical calculations. For light DM around 10 GeV, g* = 61.75, we find Xf ~ 22. 
The relic abundance is approximately given by 



6 Results 

We perform x 2 analysis to obtain upper bounds on the DM isospin violating couplings 
to light quarks. Since within uncertainties the current BESS-Polar II and PAMELA 
data can be consistent with the secondary antiproton background, we consider one-side 
exclusion limit in which only the data points below the theoretical prediction are taken 
into account, namely we calculate the quantity x 2 — Si(^i h — ^T P ) 2 / (T h ( ^ ®\ h > ^f^)' 
where $^ xp are the measured fluxes with uncertainty o"j and $* h are the theoretical 
predictions at kinetic energy T i; and set upper bounds at 95% CL. The required x 2 values 
for different degrees of freedom (d.o.f) are calculated from the standard ^-distribution. 
For instance, the required x 2 values are 64, 49.8, and 31.4 for d.o.f=47, 35, and 20 
respectively. The upper bounds obtained in this manner is more conservative than that 
using the whole data set. The constraints from the relic density are also set at 95% CL. 
For the sake of simplicity, we consider one operator at a time and ignore interference 
between the operators. 

For the operators contribute to s-wave annihilation which leads to velocity-independent 
annihilation cross sections, the relevant DM couplings to quarks are found to be tightly 
constrained by both the cosmic antiproton flux and the thermal relic density. In Fig. [3] 
the constraints on the coefficients a^ q for operator 5? = Xl^XQl^Q are shown in the 
(a5u, CLbd/cLbu) plane. The mass of the Dirac DM particle is fixed at m x = 8 GeV. The 
results show that the DAMA- and CoGeNT-favored regions are in tension with both the 
cosmic antiproton flux and the thermal relic density. 

At a 5d /a 5u = -0.89 which corresponds to f n /f p = -0.70, the DAMA and CoGeNT 
favored value is a^ u ~ 1.6 x 10 -5 GeV -2 corresponding to an annihilation cross section of 
(av ie i) ~ 3.4 x 10~ 25 cm 2 . As it can be seen in Fig. HJ for such a large cross section, the 
predicted antiproton flux is much higher than the current BESS-Polar II and PAMELA 
data and results in a huge x 2 /d.o.f=1.3 x 10 6 /35 in the "Global-Fit" model and an even 
larger one in the "Conventional" model. The upper bound set by the antiproton data at 
95% CL is a 5u < 1.7 x 10~ 6 GeV^ 2 in the "Global-Fit " model at a 5d /a 5u = -0.89, which 
is about an order of magnitude lower and corresponds to an annihilation cross section 



tth 2 



1.07 x 10 9 GeV" 1 x 



(29) 



M pl 



yfgla + 3bxf 
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Figure 3: Upper bounds on the coefficient function of a^ja^ u at 95% CL from 

cosmic antiproton flux and DM relic density. The mass of DM particle is fixed at 8 GeV. 
The favored regions and exclusion contours from various experiments such as DAMA [2], 
GoGeNT [5j, XENON [9JQ2], CDMS [HE] and SIMPLE pi] are also shown. 

(crv Te i) ~ 3.7 x 10 _27 cm 2 . The predicted antiproton fluxes are also shown in Fig. HJ In 
the two propagation models, the upper bound from the "Conventional" model is slightly 
stronger than that from the "Global-Fit" model. In Fig. HI we also plot the antiproton 
background of "Plain" model which is already higher than the data. Thus the upper 
bounds from this propagation model are expected to be much stronger. Since we are 
interested in conservative upper bounds, we do not further investigate the constraints in 
this model. In Fig. HI the upper bound from the relic density is a^ u < 6.0 x 10 -6 GeV~ 2 
at a^d/o,5u — —0.89, which is weaker than that from the antiproton flux but still in 
tension with the DAMA- and CoGeNT-favored value. 

In Fig. [5] we show the constraints on the coefficients au q for operator On = 
2m0(0*0)(gg) in (an u , a 11( i/a llu ) plane for the complex scalar DM mass fixed at = 8 
GeV. For this operator, the DAMA- and CoGeNT-favored regions are still in tension 
with the antiproton flux. At and/an u = —0.93 which corresponds to f n / f P = —0.70, 
the DAMA- and CoGeNT-favored values is an u ~ 7.9 x 10~ 6 GeV 2 , corresponding to an 
annihilation cross section of (av re \) ~ 8.5 x 10~ 26 cm 2 . The predicted antiproton fluxes 
are shown in Fig. [6] which is again much higher than the current data. The upper bound 
set by the antiproton data is an„ < 1.7 x 10~ 6 GeV 2 at a lld /a llu = -0.93 . Compared 
with the case of Os q , the constraints on the coefficients of Ow q are weaker, which is due 
to the fact that for the hadronic matrix element of scalar operator qq the Bi q factors are 
larger than that for vector operator g7 M g, which allows smaller for the same value of 
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Figure 4: Left) predictions of cosmic antiproton spectra from DM annihilation induced 
by operator 5q in the "Global-Fit" propagation model. Two cases are considered: (A) 
For a 5u = 1.6 x 10~ 5 GeV~ 2 which is favored by the DAMA and CoGeNT experiments. 
(B) For a^ u = 1.7 x 10~ 6 GeV~ 2 which is the maximal value allowed by the cosmic 
antiproton data at 95% CL. The ratio a^d/^u is fixed at —0.89 corresponding to f n / f p = 
—0.70 and the mass of DM particle is fixed at 8 GeV. The data of BESS-Polar II [31] and 
PAMELA [35] are also shown; right ) Same as left), but for the "Conventional" model. 

f P}U and results in smaller annihilation cross sections. The constraint from the thermal 
relic density is in tension with the DAMA- and CoGeNT-favored values. 

In Fig. [7J we show the favored regions and exclusion contours for operator 0^ q and 
Ow q for the coefficients corresponding to f n / f p = —0.70 in (<J p ,mDM) plane where m DM 
is the mass of the Dirac or complex scalar DM particle in the range from 5 GeV to 15 
GeV. For both operators, the DAMA- and CoGeNT-favored regions are far above the 
95% CL bounds from the antiproton flux for all the values of DM mass in this range. 
The favored regions are also in tension with the thermal relic density, especially for 0^ q . 

The operators 0\ q = XXQQ an d Ci3 g = (0* 0)g7 M g contribute to p-wave annihilation 
with cross section proportional to w 2 el . Thus they do contribute very little to the cosmic 
antiproton flux. However, their contributions to the thermal relic density cannot be 
neglected, as at freeze out the relative velocity is finite. In Fig. |8l we show the constraints 
from relic density on the coefficients a\ q and a^ q . For the operator 0\^ q one can see 
some tension between bounds set by the relic density and regions favored by DAMA and 
CoGeNT data. The constraint is not as stringent as that from the latest XENON100 
data. For the operator Oi q , the constraints from relic density is rather weak. The 
difference is again due to the different Bi q factors for these two type of interactions. 

It is straight forward to extend the discussions to Majorana fermions and real scalars. 
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Figure 5: Same as Fig. [31 but for coefficients a Uq of operator C n = 2m ^((j)* : <j))(qq). 

For the particles being its own antiparticles the primary source of the anitproton will be 
enhanced by a factor of 2 in Eq. (1241) . which may lead to more stringent constraints. 

7 Conclusions 

In summary, we have investigated the allowed values of DM-nucleon couplings in the 
scenario of IVDM for various target nuclei used in DM direct detections. We find that 
the recently updated XENON100 result excludes the main part of the overlapping signal 
region between DAMA and CoGeNT. The mass of the IVDM is restricted to a narrow 
range near 7.5-8 GeV. We have shown that whereas the effect of isospin violating scat- 
tering can relax the tensions between the data of DAMA, CoGeNT and XENON, the 
possible disagreement between some group of experiments such as that between DAMA 
and SIMPLE are not likely to be affected for any value of f n j f p . In an effective operator 
approach, we have investigated the conservative constraints on the couplings between the 
IVDM and the SM light quarks from the recent cosmic ray antiproton data and that from 
the thermal relic density. Among the four operators relevant to IVDM Oi q , 5q , Ou q , 
Oi3 q , the operators 5q and Ow q are found to be tightly constrained by the antiproton 
data and Q\% q is constrained by the relic density. Only the operator 0\ q can survive both 
the constraints while contribute to large enough isospin violating interaction required by 
the current data of DAMA, CoGeNT and XENON. The scenario of IVDM can be less 
constrained if the DM particles scatter off nuclei through some mediator particles such 
that the interaction cannot be described by effective operators. 
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Figure 6: Left) predictions of cosmic antiproton spectra from DM annihilation induced 
by operator 0\\ q in the "Global Fit" propagation model. Two cases are considered: 
(A) For a llu = 7.9 x KT 6 GeV" 2 which is favored by the DAMA and CoGeNT exper- 
iments. (B) For an u = 1.7 x 10~ 6 GeV -2 which is the maximal value allowed by the 
cosmic antiproton data at 95% CL. The ratio and/diiu is fixed at —0.93 corresponding 
to f n / fp — —0.70 and the mass of DM particle is 8 GeV. The data of BESS-Polar II [M] 
and PAMELA [35] are also shown; right ) same as left), but for the "Conventional" 
model. 
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Figure 7: The favored regions and constraints in the (cr p , m^Af) plane for operators 5q 
and 0\\ q with the coefficients corresponding to f n / f p = —0.7. 
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